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Alternative approach for description of the non-equilibrium phenomena arising in solids at a 
severe external loading is analyzed. The approach is based on the new form of kinetic equations 
in terms of the internal and modified free energy. It is illustrated by a model example of a solid 
with vacancies, for which there is a complete statistical ground. The approach is applied to the 
description of important practical problem the formation of fine-grained structure of metals during 
their treatment by methods of severe plastic deformation. In the framework of two-level two- 
mode effective internal energy potential model the strengthening curves unified for the whole of 
deformation range and containing the Hall-Petch and linear strengthening sections are calculated. 

PACS numbers: 05.70.Ln; 61.72.Bb: 61.72.Cc, 



I. INTRODUCTION 

An idea to use an additional variable for account of 
internal microstructure changes in polyatomic gases was 
proposed by Herzfeld and Rice in 1928 This idea has 
got development in other numerous researches of that 
time [2!-l5j. In 1937 Landau applied it to describe phase 
transitions in solids [6;] . In the post-war years Landau in 
cooperation with Khalatnikov and Ginzburg for descrip- 
tion of phase transitions offered some kinetic equations, 
describing the evolution of order parameter [3 [!]• In 
modern science this direction is presented by the theory 
of the phase fields [9l-[T6j. Unlike the theory of Landau 
the parameter of order in the theory of the phase fields 
is not so strict, however, this does not hinder in getting 
results which rather well coincide with the behavior of 
the real systems. 

To complete the picture, note that another direction 
in the description of such phenomena was offered in 1967 
by Coleman and Gurtin, who complemented the primary 
idea by the elements of rational mechanics fl3| ■ This di- 
rection was continued in Ref. [18,-^J. In Ref. [l^l an 
additional variable (order parameter), meaning the scalar 
density of embryonic microcracks, was introduced to de- 
scribe the initial stage of destruction of a quasi-brittle 
solid. Later in the work by Peter Van [2lj the applica- 
tion of internal state variables approach to the problem of 
destruction of cracking materials was examined. For the 
description of crackness level an additional vector vari- 
able a is introduced, which is related to thermodynamic 
conjugated variable A. 

Between both approaches there is no an impenetrable 
border, they use one and the same fundamental idea. The 
difference is that the internal state variables approach 
one comes from the principles and postulates of mechan- 
ics of continuous media, while the approach by Landau is 
based on more general energetic principles, which, how- 



' llsmet@fti.dn 



ever, need further development. Interest which is some- 
times shown in comparison of these two approaches is 
therefore clear [22]. 

A further approach of mesoscopic nonequilibrium ther- 
modynamics is currently developed for investigation of 
soft-matter problem with mesoscopic structural elements 
j25l - l28l | . The method addresses the systems entropy pro- 
duction, which is a detectable quantity at the mesoscopic 
level of condensed-matter organization; it also reveals ef- 
fectively the kinetic-thermodynamics properties and pe- 
culiarities of the corresponding dynamic matter arrange- 
ments. 

The present article is devoted to the development of 
the Landau approach, which in the present version ap- 
pears under the name of non-equilibrium evolutional 
thermodynamics (NEET). In part II, the basic postulates 
of NEET are given, which are grounded by the results at- 
tained currently. In part III the approach is applied to 
the description of defect kinetics during metal treatment 
by severe plastic deformation (SPD). In part IV summa- 
rizing conclusions are given. 



II. POSTULATES OF NEET 

The first important difference of NEET from the classic 
scheme treating the problem is the introduction of the set 
of non-equilibrium thermodynamics potentials 29, 3C|. 
Up to the present, basic, if not unique, stress was laid on 
the use of the free energy. The free energy is still a major 
function, which in energetic expression and logarithmic 
scale is a value inverse to the probability distribution 
function (PDF) of the system states |31|. A postulate is 
well known in statistical physics that the most probable 
states near a maximum of PDF, that is, automatically, 
near a minimum of the free energy [s^ are realized only. 

In this context the use of the free energy is very com- 
fortable, because the most probability (stationary) state 
can be found from the universal principle of a minimum of 
the free energy. At the same time, in structural physics 
not for all of the physical systems the PDF is certain. 
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and, consequently, the free energy on its basis. Mainly, 
one uses a phenomenological generalization of the free en- 
ergy, when the application of principle of minimum is not 
grounded statistically. Applicability of this principle is 
well proven most strictly for solids with vacancies, at the 
same time, for solids with other types of defects (disloca- 
tions, grain boundaries etc) PDF is unknown, and, con- 
sequently, applicability of this principle is not grounded 
strictly. Therefore a solid with vacancies can be a good 
model for approbation of ideas, having general character 

If the equilibrium (stationary) state of the system is 
known, in the case of small deviation of the system from 
that states, the evolution (relaxation) of the system can 
be written down as a condition that the system tends 
to the equilibrium state with a speed the greater the 
stronger the deviation. In terms of the free energy this 
condition results in the well known Landau-Khalatnikov 
evolution equation (Tj or in the case of the distributed 
systems in the Ginzburg-Landau equation Q. At the 
same time, there is no obstacle to write down the evo- 
lution of the system in terms of other thermodynamic 
potentials, for example, in terms of the internal energy. 

The internal energy is the clearest physically deter- 
mined energy of the system; it is basic in both thermo- 
dynamics and physics in general. This energy is included 
in formulation of the first law of thermodynamics; it is 
universal for both the equilibrium and non-equilibrium 
states. In addition, the generalized thermodynamic force 
is determined through it not phenomenologically, as, for 
example, for the free energy, but fully strictly within the 
framework of statistical consideration. Indeed, PDF for 
a solid with vacancies looks like (Sll . 



f{n) = CWeM^ 



U{n) 



Nlnl 



(1) 

where C is a normalizing constant, W is the thermody- 
namic probability, U{n) is the internal energy, N is the 
number of atoms in a solid, n is the number of vacancies, 
ks is the Boltzmann's constant, T is a temperature. A 
pre-exponential multiplier describes combinational, that 
is entropic, part of the distribution function, related to 
degeneration of macrostates. The exponent describes a 
restrictive part of the distribution function, related to 
overcoming the potential barriers between microstates. 
The most probable state is determined by a condition 
df{n)/dn, from where, in accordance with (1), most nat- 
urally appears a variable 



u = 



dU 

dn 



(2) 



It has a sense of the average energy of a defect (here a 
vacancy), or of the chemical potential of defects, on the 
other side, Eq. ([2|) is a typical determination of the gen- 
eralized thermodynamic force. Here it is not, however, 
postulated, but logically follows from Eq. ([T]), as a part 
of the process of determination of the equilibrium state. 



and, consequently, it must enter relaxation equation, as 
the equation describing system tendency to the equilib- 
rium state. Now the evolutional equation in terms of the 
internal energy can be written down as 
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where 7„ is a kinetic coefficient, Ue is a value of the va- 
cancy energy in the equilibrium state. 

To define the equilibrium value of the vacancy energy 
and the equilibrium density of vacancies, we all the same 
must address the condition of a maximum of PDF or a 
minimum of the free energy, and, it would seem, it is 
simpler to describe relaxation in terms of the free en- 
ergy, immediately using the Landau-Khalatnikov equa- 
tion. However, Eq. ^ is also true, and can also be 
used. Especially it can be valuable at generalizations to 
other types of defects, for which phenomenological gen- 
eralization of the free energy is problematic, and at times 
speculative, whereas the internal energy is of the univer- 
sal nature. 

The substantial difference of treating the problem in 
terms of the internal energy is that the generalized force 
^ is not equal to zero in the equilibrium state, and, 
consequently, for the internal energy the extreme prin- 
ciple does not work. This pitfall can be compensated 
as the generalized force ^ is determined uniformly for 
both the equilibrium and non-equilibrium states, if the 
equilibrium energy of the defect is known for us from 
some sources for the defined values of external control 
parameters, then it is possible to find its value for other 
parameters easily. This advantage is for the first time 
exposed in the present article. 

The second postulate, which distinguishes this ap- 
proach from traditional one, is the use of the density 
of defects as an independent thermodynamic variable in- 
stead of the configuration entropy. Note that this vari- 
able is used not in parallel with the configuration en- 
tropy, but instead of it. Mutually identical dependence 
between the density of defects and the configuration en- 
tropy can serve as a foundation for this purpose. In 
the case of a solid with vacancies this one-to-one depen- 
dence follows from the fundamental Boltzmann relation- 
ship Sc = kshiW and from definition of W in accor- 
dance with Eq. ([ij. For other types of defects this rela- 
tion is unknown, but it still must be mutually identical. 
It allows to g eneralize the 1st law of thermodynamics in 
the form 13311 



dU 



Vdijdeij 



E 

1=1 



TdS + fdS+y ifiidhi, (4) 



Eij are stress 



where V is the volume of the system, cry , 
and elastic deformations tensors, S is thermal or, for a 
solid, oscillation (not configuration) entropy, T, S are 
non-equilibrium temperature and entropy, which char- 
acterize the dynamic transitional phenomena during the 
generation and motion of the structural defects [33j . ipi, 
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Hi are energy and density of Z-kind defects, Ndef is the 
number of types of the defects. Thus the internal energy 
is a function of such independent variables as e^j, S, S 
and Hi, that is, U — U{eij,S,S,Hi). 

The first two terms in Eq. ([T]) are changes in the inter- 
nal energy due to contribution of the elastic stress field 
and equilibrium thermo-motion, the third term charac- 
terizes a part of the internal energy, arising due to non- 
equilibrium transient processes (the necessity of its ac- 
count is grounded in Ref . [s^ ) , and the last term presents 
a part of the internal energy concentrated in defect sub- 
systems. 

Relation ([5]) in the case of arbitrary number of defect 
types can be generalized as 



dU 

and evolution Eq. ([3]) accordingly 



m 

dt 



dU 

m 



fie 



(5) 



(6) 



Carrying out transformation of the Legendre type, but 
with respect to the pair of thermodynamic variables of 
ipi and Hi, we pass to a new thermodynamics potential 

Fi = U-ipiHi, (7) 

It is not hard to show that for this function a relation 

dU 



Hi = - 



d(pi ' 



(8) 



is just, that is, the new function Fi is related to the inter- 
nal energies of U in the same way as the classic (thermal 
or oscillation) free energy [37| . if the entropy is formally 
taken instead the number of defects Hi , and the temper- 
ature instead of the energy of defect (pi. But with the 
second postulate of NEET such the accordance is recog- 
nized, that is, the function Fi can be interpreted, as a 
modified free energy, when the number of defects of Hi 
is taken as an independent thermodynamic variable. 

If transformation ([7]) is done for all of types of de- 
fects, then such modified free energy will be the func- 
tion of independent variables Sij, S, S and ipi, that is, 
Fi — Fi{eij,S,S,(pi). Thus, an own argument for the 
internal energy is the defect densities, and for the modi- 
fied free energy the energy of defects, similar to the case 
of classic thermodynamics, where an own argument for 
internal energy is the entropy, and for free energy is the 
temperature. 

Classical configurational free energy Fc = U — TSc 
is always treated, as one-to-one function of the density 
of defects, while in obedience to the method of defini- 
tion (Legendre transformation) it must be the function 
of temperature. It follows that it is not a thermodynam- 
ics potential in a strict sense, and it is only useful as 
energy refiection of PDF. While the internal energy and 
the modified free energy, satisfying Eqs. JS]), ^ and ([5]), 
are real thermodynamics potentials, though dissatisfy ex- 
treme principle. 



III. TWO-LEVEL TWO-MODE MODEL OF SPD 

Now setting the dependence of the internal energy on 
its arguments, we fully determine our problem in a ther- 
modynamic sense. Let us apply the above approach for 
solution of a special problem. Producing a fine-grained 
structure of metals by severe plastic deformation (SPD) 
is presently urgent. At the initial stage of SPD there 
goes intensive generation of dislocations, then next the 
dislocations serve as a building material for the growth of 
grain boundaries, that results in a finer grain structure. 
Thus, in the processes of SPD these two types of defects 
take the main part and predetermine the two-level char- 
acter of the problem [38| . 



A. Evolution equations 

Let us consider a homogeneous problem, setting the 
internal energy as a polynomial dependence 

u = uo+ ^ ifoihi~ipiihf + (p2ihf-ip3ihf) + (pgDhghD, 

l=g,D 

(9) 

where uq, (fkh fgO are some coefficients, depending on 
the equilibrium variables of s and Eij , as control param- 
eters 

«o = ^A(4)2^M(4)' + /^«'' 
foi = foi + gi^u + IMsu? + m(4)' - (10) 

fii = fii - 2ee^j. 

For the sake of convenience we passed from the numbers 
of defects of Hi to the densities of corresponding variables 
of hi, and similarly S s, U u. 

The fourth-degree polynomial in parentheses can have 
at the positive values of coefficients fkm two maxima 
(two modes). The mode which corresponds to the lower 
value of defectiveness, in the case of dislocations I — D, 
can describe the accidental (homogeneous) distribution 
of dislocations. The mode, which corresponds to the 
higher value of defectiveness, describes dislocations be- 
longing to the cell structure in this case. We examine 
only the simplified case of the homogeneous distribution 
of dislocations, that is, ip'^o = and Lp^D = 0. 

The coefficient of (^og can be considered as a general 
surface density of the energy of regular (infinity) GB. 
From data of S.A. Firstov, for cold-roll treatment this 
energy can equal the doubled energy of a free bound- 
ary for same material. For example, it can make ap- 
proximately 2 X 2J/vn? for copper. The first term in 
this context is the own energy of the boundary without 
a contribution from other factors. It is considered, as 
a well fuU-relaxated grain boundary, that is, as a mini- 
mum possible surface energy of GB. According to A.S. 
Firstov, this energy makes 0.15 0.20 of the energy of 



4 



free-surface of the same material. That is, for copper this 
energy approximately equals 0.2 x 2 = OAJ/m^. 

Contribution of the second term ggSfj to the field of 
compressing hydrostatic stress results in the decreas- 
ing of GB energy. It is of great physical importance 
meaning that the grain boundaries are sites of density 
lack distributed along some surface. Exactly these sites 
give the highest contribution to the energy of bound- 
aries. When, due to external pressure, the volume of un- 
densitied sites diminishes, the energy of boundaries di- 
minishes, as well as the potential barriers between the 
stable neighbor states, that results in growing mobil- 
ity of the grain boundaries. If elastic deformation is 
0.002 that for copper corresponds to the level of tension 
^ ISOMPa, the constant gg is to be taken within the 
limits of 12 J/m2. For such value of constant gg the en- 
ergy of grain boundaries will diminish within the limits 
of 10 percents of its value for a relaxed boundary. 

The terms in Eq. (|10p proportional to and fig must 
give at the same level of elastic stress of ISOMPa, such 
contribution, when general effective energy of GB might 
not exceed the double energy of the free boundary (for 
copper 4J/to^). It gives conditions for choosing constants 
in the limits Xg = 0.25 ■ lOV/m^ and p,g = 0.6 • 10^ J/m^. 
The values of these constants are different because the 
effects of shear are of greater significance for structural 
rearrangement of the solid. 

Other constants are chosen issuing from the reasoning 
that the equilibrium (stationary) values of the density of 
grain boundaries were in the interval observed in exper- 
iments. For the grain boundaries these are two steady 
states with the density in a region /i^'^ = lOmm"""^ and 
^st2 _ lo^TT),-!^ where the average grain size is lOO/im 
and lOOnm, accordingly. 

The same reasoning can be repeated for dislocations. 
The minimum excess energy of dislocations, in the ab- 
sence of other factors for copper, equals approximately 
tp*^ = 5 • 10-9 J/m [39]. 

Using Eq. ([5]), which is true for both the equilibrium 
and non-equilibrium cases, one gets evolution Eqs ([6]) in 
the form 
dh]j 

= -IhoiV^loiho - llDe) + ^goihg - hge)], 



dt 



-g^ = -lh,YPgD{hD 



hoe) + Hhg - hg^)],{ll) 



where 



V = Vlg- V2g{hg + hge) + (P3g{hl + hghge + h^J. (12) 

As seen, the evolution equations do not directly de- 
pend on the parameters (pon and (pog, but can depend 
on them through the equilibrium values hz)e and hge- 
For determination of this dependence it is necessary to 
find position of the maxima of PDF. As for basic defects, 
participating in SPD, this function is unknown; it creates 
some difficulties in application of the theory for calcula- 
tion of specific systems. For overcoming the difficulties 
let us consider the method of effective potential in terms 
of the internal energy. 



B. Method of effective potential of the internal 
energy 

Let us suppose that the equilibrium energy of defect 
le weakly depends on a current value of the density of 
defects, and it can be brought under the sign of differen- 
tiation in Eq. (j6|). Then one can introduce the effective 
internal energy 



1=1 



The evolution Eq. ^ assumes a form 
dhi du 



(13) 



(14) 



Here the plus sign is selected in case if an equilibrium 
value Lpie is in the region of convexity of the internal en- 
ergy u, the minus sign is selected in the region of its 
concavity |33| . In the first case a stationary solution cor- 
responds to a maximum of the effective energy, in the 
second case a stationary solution corresponds to its min- 
imum. Formally, Eq. (jlOp realizes an extreme principle, 
as its stationary points coincide with a maximum or a 
minimum of the effective potential of the internal energy 
u. 

We take the effective energy in the same form ^ as 
the initial internal energy u with the same coefficients of 
presentation ([TUl) with the only difference that the equi- 
librium energy of ^pie is included in coefficient of ipoi , that 
is, ifQi > (poi — (file and u > u. Then, the set of evolution 
Eqs (jlip can be written in the explicit form 



dhE 
dt 



^^9 I I, 

-Q^ = -IhS'P^g - Vlghg 



'fi2ghl - fSghl + fgDho), 



Results got directly from the solution of the set of evo- 
lution Eqs (ITT|) . and by the method of effective poten- 
tial of the internal energy ([15]), coincide, if relations be- 
tween the coefficients of internal energy expansion hge — 

{ipOg - V'ge)/Vlg, Vlg » V2ghge and Lfig » ip3g{hge)'^ 

are fulfilled. 

Coming from the above analysis, such set of param- 
eters and coefficients was accepted for calculations A — 
fi = 2.08-10i°Pa, ip*„jj = 5-10-Vto-i, ip^jj = 10 
2 • 10-^Jm-\ an = 3.3 ■ 10-4 Jm-i, 



^19 

3 ■ 10~20Jto, gg 



SD = 

f_2g = 
12Jm-2 X ^ 



9D 

6 • 10-23 Jto, ip*g = 0.4Jm-2, = 3 • 10"^ Jm^^ 
5.6 • 10-13 J, <^3g 

2.5 • lO^Jm-^, fig = 6- 10Vto-2, Cg = 3.6 • 10^'^Jmr\ 

P>gD = 10-16 J. 

Other coefhcients in expressions for the internal en- 
ergy are considered to be zero. A time step in numerical 
calculations is r = 0.67 • 10-^, kinetic coefRcients are 
Iho = 5 • IQ^^Jms, 7/i — 10^ Jm-^s. The calculation of 
system evolution with these parameters and coefficients 
is shown in fig. [T] 
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FIG. 1. Regularities of defect-formation during SPD: a) ki- 
netics of defects: 1 - density of grain boundaries; 2 - density 
of dislocations; b) unified curve of strengthening: 1 - region 
of the Hall-Petch law; 2 - region of the linear law of strength- 
ening 



From fig. [T]it is evident, that kinetics of grain bound- 
aries and dislocations during the structural phase tran- 
sition is closely correlated. At the first stage a grovifth 
in the number of dislocations initiates growth of grain 
boundaries and provokes the beginning of the structural 
phase transition. During the structural phase transition, 
when the density of dislocations has already gone on a 
stationary plateau, vice versa, the growth of grain bound- 
aries provokes the growth of the density of dislocations. 
Thus, in this area, dislocations follow the grain bound- 
aries in repeating the shape of curve of the structural 
phase transition, but to more weak extent. 



C. Strengthening curves 

As known the law of strengthening results from disloca- 
tion mobility decrease due to braking by different defects, 
as well as by dislocations from other slide planes. At the 
dislocation level the law of strengthening is described by 



Taylor relation [13, Ei: 

r = afib\/ ho, (16) 

where r is the shear stress; a is a coefficient, which takes 
on a value from the interval 0,2 — 1,0; /i is the shear 
modulus; b is the Burgers vector; ho is the density of 
dislocations (Here the modified denotation is used, which 
is general in a multilevel system of defects). 

If GB are formed directly due to an outcoming of dislo- 
cations, in this case one can be limited by relations (|16l) . 
and the Hall-Petch law at the grain level is got, as simple 
consequence from this relation acting at the dislocation 
level. 

In the theory of NEET deformation is a control pa- 
rameter which in the case of shear deformation relates to 
the stress by a simple dependence r = /le*^ so, in terms 
of the theory the law of strengthening looks like 

= ab^/h^. (17) 

This dependence can be used in NEET, as an addi- 
tional relation to energetic and kinetic relations written 
above. 

In the case of GB, a mode, which corresponds to a 
lower value of defectiveness, describes a coarse-grained 
structure, a mode, which corresponds to the greater value 
of defectiveness, describes a fine-grained structure. A 
possibility of forming different modes of the same defect 
is related to the microscopic mechanisms of deformation. 
In the case of grains they can be related to the circum- 
stance that at the initial stage an increase of the general 
surface of GB is the effective mechanism of energy dis- 
sipation. The contribution of triple junctions can be ig- 
nored at this stage. The situation changes substantially, 
when the average size of grains decreases to 100 nm. In 
this case, triple junctions can give a considerable contri- 
bution to the energy of boundaries, which can result in 
the formation of a new maximum in this region. 

Triple junctions can be considered, as a specific type 
of defect, but as it is topologically attached to the grain 
boundaries, they can be integrated and considered, as 
one defect with a somewhat more difficult dependence of 
its energy on the number of defects. 

The deformation stages during SPD are demonstrated 
by changes in the character of strengthening law at dif- 
ferent stages of the process. At the first stage (area 1, fig. 
[T]d) the law of strengthening can be approximated by the 
Hall-Petch law, if expressed through elastic deformations 

£^ = e-^o - A^g. (18) 

where the constants are e^^ = —0.075 and A = 0.119mm. 
The negative sign is taken because at compression that 
is in conditions typical of SPD the elastic deformation 
is negative. As seen from the picture, the elastic defor- 
mation changes within the limits of 0.1 0.5 procents, 
that in view of the value of shear modulus is in region of 
values for real materials. 
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In region 2, which corresponds to the most rapid phase 
of deformation during SPD, the law of strengthening can 
be approximated by a Unear dependence 

e"" = e"^ - Bhg. (19) 

with constants e^^ — —0.27 and B = O.OOSOmm^. Just 
the same character of strengthening law change depend- 
ing on deformation stage was noted in Ref. [40]. 



IV. FINAL REMARKS 

In the article, an alternative approach of non- 
equilibrium evolutional thermodynamics is considered 
and all features of the approach are demonstrated by 
the model example of solids with vacancies. We derive 
the system of kinetic equations in terms of the internal 
energy as the most fundamental thermodynamic poten- 
tial. The generalized thermodynamic force immediately 
follows from differentiation of the probability distribu- 
tion function during the procedure of finding the most 
probable state. In the equilibrium (or stationary) state 
such force is not zero, and associated with the equilib- 
rium energy of defect in the equilibrium state. 

The second feature which distinguishes this approach 
from traditional one is in using the density of defects 
as an independent thermodynamic variable instead of 
the configuration entropy. This allowed to modernize 
the writing of 1st law of thermodynamics, by adding 
to it an entropy terms in the form of product of the 
defect energy and the increment of the defect density 
(the last term in Eq. In addition the first law 

of thermodynamics in this approach is written with dy- 
namical transient phenomena taken into account in the 
form thermodynamic processes, to this end, the concepts 
of non-equilibrium temperature and non-equilibrium en- 
tropy were introduced (term before the last one in Eq. 
(HJ). Such a formulation can not be treated as funda- 



mental, as the temperature and the entropy introduced 
in such a way will be determined by statistics of transient 
phenomena, which is not universal but strongly depen- 
dent on character of an external infiuence. At the same 
time, it can useful as an approximated relation for the 
solution of concrete problems. 

Inclusion in the modified 1st law of thermodynamics 
of additional terms, describing internal non-equilibrium 
processes (they can be considered as analogues of the in- 
ternal variables or the order parameters introduced by 
Landau), really, extends the dimensionality of the prob- 
lem. The increase of problem dimensionality makes it 
simply certain. At the same time, we can introduce the 
generalized concept of the system state, depending on 
both equilibrium e^p S and non-equilibrium S and Hi 
variables. 

Approach of NEET is applied to simulate the extraor- 
dinarily important for practical applications problem, the 
refining of grain structure of metals by severe plastic de- 
formation methods. Within the framework of the two- 
level and two-mode approximation the system of kinetic 
equations, which includes the mutual kinetics of disloca- 
tions and grain boundaries, is deduced. For practical cal- 
culations the comfortable phenomenological model of the 
effective thermodynamic potential of the internal energy, 
which realizes the usual concept of the extreme state, is 
introduced. 

Within the framework of this model the unified curve 
of strengthening, which has stage-like character, is con- 
structed over the whole of the deformation interval. At 
the initial stage the curve of strengthening can be ap- 
proximated by square-root dependence close to the Hall- 
Petch law, and at the finishing stage, it can be easily 
approximated by linear dependence. 
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